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LUBRICATION OF RIGID ELLIPSOIDAL SaiDS 

There have been relatively few theoretical studies of the 
lubrication of rolnng contacts compared with the attention 
given to plain thrust and journal bearings. Certainly a great 
deal of development effort has been put into gear, ball, and 
roller bearing lubrication since all are machine elements of 
great technological and commercial importance. But understand- 
ing of their lubricating mechanism is still far from complete. 
The huge pressures involved in the elliptical contact generated 
between a rolling element and a track influence in a rather com- 
plicated way the properties of the lubricant and the shapes of 
the contacting surfaces. 

The first requirement therefore is to develop a basic solu- 
tion to the problem of the lubrication of rigid ellipsoidal 
solids with an isoviscous, incompressible fluido A solution to 
this problem is presented in this chapter, and the results pro- 
vide a foundation for the full appreciation of the elastohydro- 
dynamic theory of elliptical contacts to be presented later. 

The influence of the variation of viscosity with pressure and 

*Published as Chapter 6 in Ball Bearing Lubrication by 
Bernard J. Hamrock ana Duncan Dowson, John Wiley & Sons, 
Inc., Sept. 1981. 
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the effect of the deformation of surfaces under pressure are 
considered in Chapter 7. 

Most studies of the hydrodynamic lubrication of nominal 
point and Une contacts have concentrated on minimum-film- 
thickness predictions for either a ball near a plane or a cylin- 
der near a plane in which side leakage is neglected. But the 
full range of geometries between the two extremes has only re- 
cently been studied by Brewe, et al. (1979). Kapitza (1955) 
presented an early and elegant analysis of this problem in which 
he generated a minimum-film-thickness fonnula for a nonrotating 
sphere floating in a sea of lubricant above a plane surface 
sliding at a given velocity. This analytical solution repre- 
sents a rare and outstanding example of a successful mathemati- 
cal approach to the solution of the second-order differential 
equation presented by Osborne Reynolds. Solutions taking ac- 
count of side leakage were developed analytically for the spe- 
cial case of a ball on a plane by using a clever substitution. 
One disadvantage of the Kapitza approach is that it adopted the 
half-Sonimerfeld boundary condition, which violates the require- 
ment of flow continuity at the cavitation boundary. However, 
the effect of this assumption on film thickness prediction is 
not always serious. The Brewe, et al. (1979) work is used ex- 
tensively in this chapter because it uses the Reynolds cavi- 
tation boundary conditions and is applicable to the complete 
range of geometries of the contacting solids. 
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The materia! presented in this chapter not only provides an 
introduction to classical hydrodynamic lubrication as it relates 
to nonconformal contacts, but it also introduces a minimum- 
f ilm-tnickness formula that has a number of direct applica- 
tions: circular-arc thrust bearing pads; 'nidustrial coating 
processes in which paint, emulsion, or protective coatings are 
applied to sheet or film materials passing betweens rollers; and 
very lightly loaded cylindrical roller bearings. 

6.1 Equation for Pressure Distribution 

From the previous chapter the appropriate form of the 
Reynolds equation when the viscosity is assumed to be constant 
is given as 

Is) + IS 


where 


u 
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It is convenient to nondimensionalize with respect to the effec- 
tive radius R^; that is, 




and IP 


pR 

X 

T)o“ 


( 6 . 1 ) 


in terms of these dimensionless variables equation (5.45) becomes 
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The film thickness between two rigid ellipsoidal solids as 
shown in Figure 2,18 can be written as 

h - hQ + S (6.3) 

where 

S a geometrical separation of the solids 
hg « central (minimum) film thickness 

Using equation (2,35) and rewriting equation (f>.3) in dimension- 
less form give 


where 


H s. K, 



(6,4) 


“a-r (6.5) 

ho 

"o’r 

X 

In equation (6,4) the parabolic approximation is used in 
defining the geometrical separation of the undeformed solids. 

The solution of the Reynolds equation (6.2) consists of a 
homogeneous and a particular solution; thqt is, 


Ph + Pp 


(6.7) 
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for which is a solution to the homogeneous equation and 
the condition that Pj^ » -Pp at the bouidaries leads to 


-L 

9X 



( 6 .«) 


Kapit?a (1955) was the first to recognize that the particu- 
lar solution for the pressure is simply proportional to X/H^, 
or 



where 

“ 1 + 2/3a (6.10) 

a 

In the preceding equation is the side-leakage factor es- 
tablished by Archard and Cowking (1965-66) and can be verified 
by inserting Pp back into equation (6.2). If we define 
P[^(X,Y) s 4XjjQ(X,Y) by using equation (6.7), we can express the 
full solution as 

" ■ * 'i) 

In general the homogeneous solution P|^ is an unknown func- 
tion of X and Y. Consequently, the pressure distribution 
must be determined numerically for a given speed, viscosity, 
geometry, and film thickness. The numerical solution is 


normally achieved by the Gauss-Seidel Iterative method witli 
overrelaxation. 

A variable-mesh nodal structure like the one shown in Fig- 
ure 6.1 Is used to provide close spacing In and around the pres- 
sure peak. In this figure the inlet is to the left and the out- 
let is to the right, The variable mesh helps to minimize the 
errors that can occur because of large gradients in the high- 
pressure region. The grid spacing of the coordinates X and 
Y varies depending on the anticipated pressure distribution. 
That is, for a very highly peaked and localized pressure distri- 
bution, the dimensionless fine mesh spacing is normally about 
0.002, and the coarse mesh about 0.1. For a relatively flat 
pressure distribution, the fine mesh would be about 0.005, and 
the coarse mesh about 0.13. 

6.2 Boundary Conditions 

At the inlet the pressure is taken as zero at a distance 
sufficiently far from the center of the contact to give fully 
flooded conditions. Similarly, at a sufficient distance from 
the center of the contact to the sides of the nodal structure, 
the pressure is also taken to be zero. The boundary condition 
in the outlet is not as straightforward since it is necessary to 
take account of cavitation. 
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In thf, outlet region there is a tendency to form subambient 
pressures, which lead to disruption of the lubricating film by 
cavitation. The usual form of cavitation in lubricating films 
is the liberation of dissolved gases. Mineral oils contain be- 
tween 8 and 10 percent of dissolved air. When the pressure in 
the oil film falls below ambient, some of the air is liberated 
in the form of bubbles. This tends to maintain the oil film 
pressure near the level of the saturation pressure. For most 
lubrication conditions the saturation and ambient pressures will 
be almost equal. These observations suggest that the pressure 

I 

in the cavitated region of lubricating films will be anproxi- 
mately constant and near to the atmospheric or ambient pre.ssure. 

The approach adopted by Kapitza (1955) in defining the out- 
let boundary condition was to ignore the negative pressures, 
that is, to employ the half-Sommerfeld boundary condition. 
Kapitza' s solution has the appeal of simplicity, but it does not 
satisfy the flow continuity conditions at the cavitation 
boundary, namely, that the pressure gradient normal to the cavi- 
tation boundary must be zero. To insist that P = aP/aN = 0 
at X = 0 would be overspecifying the problem mathematically. 

It is possible, however, to insist that P « aP/aN = 0 at the 
cavitation boundary and to locate the interface in such a posi- 
tion that this well-known Reynolds boundary condition is 
satisfied. 


6.3 Load Capacity 


Once the pressure distribution for the appropriate cavlta-* 
tion boundary conditions has been determined numerically, we can 
express tne load capacity as 

P dx dy (6,12) 


By making use of equation (6.11) we can write this equation in 
dimensionless form as 


r 





(6.13) 


The centra] film thickness can be isolated from the integrand by 
introducing the following transformations: 


X - 


(6.14) 


If we assume the homogeneous solution to transform In the 
same manner as the particular solution, we obtain 


F “ SXjjfiQUR 



-X. 


1 + X? + yf 

t ti 


T + Q(Xt.Xt) 


dXt dY^ 


(6.15) 
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Kapitza (1955) refers to this integral as the reduced hydrody- 
namic lift L^. Thus 

dXj. dY^ ( 6 . 16 ) 


6.4 Film Thickness Formula 

The reduced hydrodynamic lift was found by Kapitza (1955) 
to equal it/2 by as^acving Q « 0 and integrating over the 
half-space of positive pressures. For the Reynolds boundary 
conditions the limits depend on the shape of the cavitated re- 
gion and hence on the geometry. Consequently we seek an addi- 
tional geometrical effect to modify and generalize Kapitza's 
solution. The central (minimum) film thickness can be expressed 
as a function of the load, spee;j, geometry, and fluid viscosity 
by rearranging equation (6.15) and writing 

/A.rinuRv 

«n,ln - »0 - “a (-4-^ 

The ratio of dimensionless speed to dimensionless load can be 
defined as 



•X,. 


(l + + A- 


+ QCX^.Yj.) 
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( 6 . 18 ) 



and equation (6.17) becomes 

2 

“a (^b^t m) (6*1^) 

The integrand of equation (6,16) thus becomes a function of the 
geometry represented by and the central film thickness 
Hq. Consequently ■ L^(ag,HQ)> and this results in a trans- 
cendental equation for Hq, 

It is necessary to determine as a function of the 
geometry alone, or 

Lj. - L^(aa) if xl « (6. 2:0) 

Once the hydrodynamic load-carrying capacity F has been 
obtained from the numerically determined pressure distribution, 
can be evaluated for various geometries. A curve fit can 
then be used to show the effect of geometry on L^. The fol- 
lowing relationship provides a good representation of the numer- 
ical results; 

Lj. « 0.131 tan“^ + 1.683 (6.21) 

Inserting equation (6.21) into equation (6.19) gives the follow- 
ing general expression for minimum film thickness between any 
rigid ellipsoidal solids, ranging from two spheres in nominal 
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point, contact to Infinitely long cylinders in nominal line 
contact: 



6.5 Comparison Between Different Theories for the Lubrication 
of Rigid Ellipsoidal Solids 


The minimum-film-thickness equation derived by Kapitza 
(1955) using the half-Sommerfeld cavitation boundary condition is 

<Wk ■ 1^8 “a (i t)' 

From equations (6.19) and (6.23) we can write 



By using equations (6.21) and (6.24) we find that the film thick- 
ness as obtained from equation (6.22) is 11 to 21 percent 
greater than that obtained from Kapitza's (1955) solution 
(equation (6.23)), with the smallest difference occurring for a 
ball-on-plane contact. The alteration of the pressure distribu- 
tion due to the Reynolds cavitation boundary conditions is re- 
sponsible for this influence of contact gometry on minimum film 
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thickness. These differences between the results of Kapitza 
(1955) and Brewe, et al. (1979) are illustrated graphically in 
Figure 6,2. For the entire range of the radius ratio R„/R„ 
the dimensionless film thickness 1s higher for the Brewe, et al. 
(1979) results than for the Kapitza (1955) results. 

Both the Kapitza (1955) analysis and the numerical solu- 
tions of Brewe, et al. (1979) resulted in an exponent of 2 
for U/W in the dimensionless film thickness equation. Dalmaz 
and Godet (1973) also used th«? *\eynolds boundary conditins for a 
ball-on-plane configuration, and they reported a comparable ex- 
ponent of 1.77. This lower exponent has been discussed by 
Bre^’ie, et al. (1979), and it appears to be due to starvation 
effects resulting from the inlet boundary condition for both the 
analytical and experimental (Dalmaz and Godet, 1973) situa- 
tions. The sophisticated apparatus used by Dalmaz and Godet 
simultaneously measured the film thickness, traction, and speed 
between a steel ball and a glass plate in pure sliding. A 30-mm- 
diameter ball, which turned around one of its axes, was immersed 
in an oil bath and thus carried the lubricant into the contact 
area through viscous lifting. An optical interferometric tech- 
nique was used to measure the oil film thickness. 

Experimental results obtained by Dalmaz and Godet (1973) 
under lightly loaded, isoviscous conditions for pure sliding of 
a ball on a plane are compared with the theoretical results ob- 
tained from equation (6.22) in Figure 6.3. The grouping of 


dlniansionless parameters in this figur^ is the same as that used 
by Dalmaz and Godet (1973) and was first introduced by Thorp and 
Gohar (1972). The theoretical results obtained from equation 
(6.22) are in excellent agreement with the experimental data for 
the lower range of Hq/WG, but at higher speeds the measured 
film thicknesses fall below the theoretical predictions, prob- 
ably as a result of lubricant starvation in the inlet region. 

For comparison, predictions based on the theory of Kapitza 
(1955) and the calculations of Dalmaz and Godet (1973) have also 
been included in Figure 6.3. 

6,6 Pressure Distribution Between Ellipsoidal Solids 

Three-dimensional contour plots of the pressure distribu- 
tion for values of of 1.00 and 36.54 are presented in 
Figure 6.4. The shape of the cavitation boundary is clearly 
evident in this figure. As becomes large, the cavitation 
boundary tends to straighten out and is accompanied by de- 
creasing changes in L^. The scale along the Y axis in Fig- 
ure 6.4 has been magnified about three times to improve the 
resolution. Consequently the differences in the shapes of the 
cavitation boundary are actually subdued in Figure 6.4. 

Isobar plots for three radius ratios of 25.29, 8.30, 
and 1.00 are shown in Figure 6.5. The center of contact is rep- 
resented by an asterisk. The pressure peak builds up in the 
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entrance region, wnicii is located to the left of the center of 
contact and 1s indicated by a cross. Since the isobars in each 
case are evenly spaced, the pressure gradients can be easily 
envisaged. Note that, as the radius ratio ct^ increases, the 
steeper pressure gradients are found predominantly in the roll- 
ing direction. This implies that the amount of side leakage 
decreases as increases. A decrease in side leakage is 
reflected in an increase in the value of X|^. For nominal 
line contact « i, and for the largest value of Og (25.29) 
recorded in Figure 6.5 the corresponding value of Xj^ is 
0.974. 

6.7 Closure 

The influence of geometry on the hydrodynamic lubrication 
of rigid, ellipsoidal soli,; has been investigated in this chap- 
ter. The study has been restricted to conjunctions fully im- 
mersed in lubricant (i.e., fully flooded). The effect of geome- 
try on film thickness was determined numerically by varying the 
radius ratio from 1 (a ball-on-plane or ball-on-ball con- 
figuration) to 36 (a ball in a conforming groove). Pressure- 
viscosity effects were not considered, It was found that the 
minimum film thickness had the same speed, viscosity, and load 
dependence as found by Kapitza in his classical analytical solu- 
tion to the problem, When the Reynolds cavitation boundary con- 
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dition was incorporated in the analysis, an additional geometri- 
cal effect was introduced into the film ^riickness equation. The 
derived film thickness equations can be compared as follows: 
Reynolds boundary condition: 



Half-Sommerfeld boundary condition: 

■ 128 “a (¥i)' 

With the Reynolds boundary condition the minimum film thickness 
has been found to be 11 to 21 percent greater than that obtained 
when using the half-Sommerf eld boundary condition adopted by 
Kapitza. 
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SYMBOLS 


A 


A^, B*, C*, 
D*. L*, M* , 


a 

a 


B 

b 

b 

C 


’ 1 ’ ■ 


’-8 


c 


c 


D 

D 

D 


De 


d 

d 




e 




d 

d 


i 

0 


constant used In equation (3.113) 
relaxation coefficients 

2 

drag area of .ball^ m 

semimajor axis of contact ellipse, ni 

a/ 2m 

total conformity of bearing 
semiminor axis of contact ellipse, m 
b/2m 

dynamic load capacity, N 
drag coefficient 
constants 

19,609 N/cm^ (28,440 Ibf/in^) 

number of equal divisions of ?emimajor axis 

distance between race curvature centers, m 

material factor 

defined by equation (5.63) 

Deborah number 
ball diameter, m 

number of divisions in semiminor axis 
overall diameter of bearing (Figure 2.13), m 
bore diameter, m 
pitch diameter, m 

pitch diameter after dynamic effects have acted on ball, m 
inner-race diameter, m 
outer-race diameter, m 
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modulus of elasticit)/, N/m^ 

effective elastic modulus, 2 

internal energy, m^/s^ 
processing factor 



* N/m^ 


p X/fc 

elliptic integral of second kind with modulus (1 - l/k‘) 

approximate elliptic integral of second kind 

dispersion exponent 

normal applied load, N 

normal applied load per unit length, N/m 

lubrication factor 


integrated normal applied load, N 
centrifugal force, N 

maximum normal applied load (at » 0), N 

applied radial load, N 

applied thrust load, N 

normal applied load at angle \l>, N 

1/2 

elliptic integral of first kind with modulus (1 - 1/k^) 

approximate elliptic integral of first kind 

race conformity ratio 

rms surface finish of ball, m 

rms surface finish of race, m 

dimensionless materials parameter, uE 

fluid shear modulus, N/m^ 

haroness factor 

gravitational constant, m/s^ 
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9e 

9V 

H 

a 



u 

min,r 

H . 
rmn,s 


H 


c 


u 

min 


u 

min 


n 

h. 




dimensionless elasticity parameter, 
dimensionless viscosity parameter, GW^/U^ 
dimensionless film thickness, h/R 

A 

dimensionless film thickness, H{W/U)^ •« F^h/u^no^x 
dimensionless central film thickness, 
dimensionless central film thickness for starved 
lubrication condition 
frictional heat, N m/s 

dimensionless minimum film thickness obtained from EHL 
elliptical-contact theory 

dimensionless minimum film thickness for a rectangular 
contact 

dimensionless minimum film thickness for starved 
lubrication condition 

dimensionless central film thickness obtained from 
least-squares fit of data 

dimensionless minimum film thickness obtained from 
least-squares fit of data 

dimensionless central-film-thickness - speed parameter, 
H,U-0-5 

dimensionless minimum-film-thickness - speed parameter, 
H 

mm 

new estimate of constant in film thickness equation 
film thickness, m 
central film thickness, m 
inlet film thickness, m 
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m 


imn 


0 

0* 

J 

K 

k 

k 

k 

kf 

L 

••a 

“■t 

L , , . . , 

^10 


•50 


a 

T 

M 


film thickness at point of maximum pressure, where 
dp/dx «• 0, m 

minimum film thickness, m 
constant, m 

diametral interference, m 

0 

ball mass moment of inertia, m N s 
integral defined by equation (3.76) 
integral defined by equation (3.75) 
function of k defined by equation (3.8) 
mechanical equivalent of heat 
polar moment of inertia, m N s^ 
load-deflection constant 
ellipticity parameter, a/b 
approximate ellipticity parameter 
thermal conductivity, N/s *C 
lubricant thermal conductivity, N/s "C 
fatigue life 
adjusted fatigue life 

reduced hydrodynamic lift, from equation (6,21) 
lengths defined iio Figure 3.11, m 
fatigue life where 90 percent of bearing population will 
endure 

fatigue life where 50 percent of bearing population will 
endure 

bearing length, m 

constant used to determine width of side-leakage region 
moment, Nm 
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M 

M 

M 


g 

p 

s 


m 


m* 


m 

m 


m, 


W 


N 

n 

n* 

n’ 

P 

Po 


Hz 


max 


IV, as 


m 


Sfjyroscopic moment, Nm 

n 75 

dimensionless load-speed parameter, WU ' 
torque required to produce spin, N m 
mass of Pali, N s^/m 

dimensionless inlet distance at boundary between fully 
flooded and starved conditions 
dimensionless inlet distance (Figures 7.1 and 9.1) 
number of divisions of semimajor or semiminor axis 
dimensionless inlet distance bounoary as obtained from 
Wedeven, et al. (1971) 
rotational speed, rpm 
number of balls 
refractive index 

constant used to determine length of outlet region 

dimensionless pressure 

dimensionless pressure difference 

diametral clearance, m 

free endplay, m 

p 

dimensionless Hertzian pressure, N/m 

p 

pressure, N/m 

p 

maximum pressure within contact, 3F/2»ab, N/m 

^ p 

isoviscous asymptotic pressure, N/m 
solution to homogeneous Reynolds equation 
thermal loading parameter 

dimensionless mass flow rate per unit width, qno/pQE'R^ 
reduced pressure parameter 

p 

volume flow rate per unit width in x direction, m /s 
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R 

Ra 

R, 


r 

r 


ax 


’ ^by J 


ay’ 

»'c’ ^c’ " 

"'s’ ®s’ 

F 

S 

S* 

So 

S 

T 

t 

■^m 

^b* 

AT* 

T, 


t 

t. 


2 

volume flow rate per unit width in y direction, nr/s 
curvature sunii, m 

arithmetical mean deviation defined in equation (4.1), m 
operational hardness of bearing material 
effective radius in x direction, m 
effective radius in y direction, m 
race curvature radius, m 

radii of curvature, m 

cylindrical polar coordinates 

spherical polar coordinates 

defined in Figure 5.4 

geometric separation, m 

geometric separation for line contact, m 

empirical constant 

shoulder height, m 

‘''O^Pmax 

tangential (traction) force, N 
temperature, *C 
ball surface temperature, "C 
average lubricant temperature, *C 
ball surface temperature rise, "C 

viscous drag force, N 
time, s 

auxiliary parameter 

velocity of ball-race contact, m/s 
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velocity of ball center, m/s 

dimensionless speed parameter, nQd/E'Rj^ 

surface velocity In direction of motion, (u^ + u^)/2, m/s 

number of stress cycles per revolution 

sliding velocity, u^ - Ujj, m/s 

surface velocity in transverse direction, m/s 

dimensionless load parameter, F/E'R^ 

surface velocity in direction of film, m/s 

dimensionless coordinate, x/R 

A 

dimensionless coordinate, y/R^ 

dimensionless grouping from equation (6.14) 

external forces, N 

constant defined by equation (3.48) 

viscosity pressure index, a dimensionless constant 

coordinate system 

pressure-viscosity coefficient of lubrication, m^/N 

radius ratio, R^/R^ 

contact angle, rad 

free or initial contact angle, rad 

iterated value of contact angle, rad 

curvature difference 

viscous dissipation, N/m^ s 

total strain rate, s“^ 

elastic strain rate, s“^ 

viscous strain rate, s“^ 
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s 

6 * 


'D 


& 


6 


c 


c 


n 

0 

A 

X 


flow angle, deg 

total elastic deformation, m 

0«"1 

lubricant viscosity temperature coefficient, C 
elastic deformation due to pressure difference, m 
radial displacement, m 
axial displacement, m 
displacement at some location x, m 
approximate elastic deformation, m 
elastic deformation cf rectangular area, m 
coefficient of determination 
strain in axial direction 
strain in transverse direction 
angle between ball rotaiiOiial axis and bearing 
centerline (Figure 3.10) 
probability of survival 
absolute viscosity at gauge pressure, N s/m*^ 
dimensionless viscosity, n/nQ 

2 

viscosity at atmospheric pressure, N s/m 
6.31x10“^ N s/m2(0.0631 cP) 
angle used to define shoulder height 
film parameter (ratio of film thickness to composite 
surface roughness) 

equals 1 for outer-race control ana 0 for inner-race 
control 

second coefficient of viscosity 
Archard-Cowking sidt»-leakage factor, (1 + 2/3 
relaxation factor 
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V 


M* 

V 

t 

P 

■p 

PQ 

0 

T 

”^0 

T 


T 


e 


$ 


0 

h 


«0 


“b 


coefficient of sliding friction 
"pR 

Poisson's ratio 

divergence of velocity vector, (tu/3x) + (8V/ay) + (aw/az), 
lubricant density, N s^/m^ 
dimensionless density, p/pq 

? 4 

density at atmospheric pressure, N s‘-/m 
normal stress, N/m^ 
stress in axial direction, N/m 
shear stress, N/m^ 

2 

maximum subsurface shear stress, N/m 
shear stress, N/m^ 

2 

equivalent stress, N/m 

2 

limiting shear stress, N/m 

ratio of depth of maximum shear stress to semiminor axis of 
contact ellipse 
Ph3/2 

auxiliary angle 
thermal reduction factor 
angular location 
limiting value of 

absolute angular velocity of inner race, rad/s 
absolute angular velocity of outer race, rad/s 
angular velocity, rad/s 
angular velocity of ball-race contact, rad/s 
angular velocity of ball about its own center, rad/s 
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Wg angular velocity of oall around shaft center, rad/s 

Wg ball spin rotational velocity, rad/s 

Subscripts: 


a 

solid a 

b 

solid b 

c 

central 

be 

ball center 

IE 

isoviscous-elastic regime 

IR 

isoviscous-rigid regime 

1 

inner race 

K 

Kapitza 

min 

minimum 

n 

iteration 

0 

outer race 

PVE 

piezoviscous-elastic regime 

PVR 

piezoviscous-rigid regime 

r 

for rectangular area 

s 

for starved conditions 

x,y ,2 

coordinate syste.ii 

Superscript: 


(~) 

approximate 
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Figure 6. 1. - Variable nodal structure used for numerical 
calculations. 
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Figure 6,2, - Film thickness results of Kapitza (195S) and Brewe, et al, (1979) for range 
of R^/Ry and U/W « lxlO'3 
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